Interplay between magnetism and superconductivity in Fe-pnictides 



00 

o 
o 

(N 
o 

Q 



O 

o 
I 

a.: 



O 

o 



00 

o 



X 



A. B. Vorontsov, M. G. Vavilov, and A. V. Chubukov 
Department of Physics, University of Wisconsin, Madison, Wisconsin 53706, USA 

(Dated: December 12, 2008) 

We consider phase transitions and potential co-existence of spin-density-wave (SDW) magnetic 
order and extended s-wave (s""") superconducting order within a two-band itinerant model of iron 
pnictides, in which SDW magnetism and superconductivity are competing orders. We show that 
depending on parameters, the transition between these two states is either first order, or involves 
an intermediate phase in which the two orders co-exist. We demonstrate that such co-existence is 
possible when SDW order is incommensurate. 
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Introduction. Iron-based pnictide superconductors - 
oxygen containing 1111 materials RFeAsO (R = La, Nd, 
Sm) and oxygen free 122 materials AFe2As2 (A = Ba, 
Sr, Ca) are at the center of experimental and theoret- 
ical activities at the moment because of high potential 
for applications and for the discovery of new mechanisms 
of superconductivity. Most of parent compounds of Fe- 
pnictides are magnetically ordered. Upon doping, mag- 
netism disappears and superconductivity emerges, but 
the nature of this transition remains unclear. Some ex- 
periments on fluoride-doped 1111 materials indicate that 
the transition is first-orderfi some show behavior more 
consistent with a quantum-critical point separating mag- 
netic and superconducting (SC) states^., and some show 
a co-existence of magnetism and SC in both classes of 
materials i^i^i^ 

It is by now rather firmly established that Fe-pnictides 
are metals in a paramagnetic phase for all dopings, with 
two sets of (almost) doubly degenerate Fermi surface 
(FS) pockets - a hole pocket centered at (0,0) and an 
electron pocket centered at tt = (tt, tt) in the folded Bril- 
louin zone. To a good approximation, hole and electron 
FS are circular and at zero doping have nearly identical 

c^^oc 6.7.8.9.10.11 



Like in chromium (Ref. Il2i). this nesting 
geometry is favorable to a spin-density-wave (SDW) or- 
dering at TT as the corresponding susceptibility logarith- 
mically diverges at T = and a small repulsive inter- 
action in the particle-hole channel at momentum transfer 
TT already gives rise to an SDW instability at T = T^. If 
the interaction is attractive in a SC pairing channel, then 
the SC pairing susceptibility also diverges logarithmically 
at T = and the system becomes a SC at T = Tc, unless 
magnetism interferes. 

Previous studies of an itinerant model of small elec- 
tron and hole FS have found that the same interac- 
tion, inter-band Josephson-typc pair hopping, gives rise 
to an SDW order and to superconductivity with ex- 
tended s-wave {s~^) symmetry of the SC order parameter 
(A(A;) oc coskx + cos ky in the folded Brillouin zone)r^^ 
leading to competition between the two orders. The full 
interactions in SDW and s"*~ channels also involve inter- 
band forward scattering and intra-band Hubbard inter- 
action, respectively, and the full interaction is larger in 
the SDW channelji^ Then at zero doping, which we as- 



sociate with near-perfect nesting, the highest instability 
temperature is that of an SDW state. At a nonzero dop- 
ing X, nesting is destroyed (either hole or electron pocket 
gets relatively larger), and SDW susceptibility no longer 
diverges. Magnetic Ts{x) then goes down with doping 
and above a particular value of x, the first instability 
upon cooling is into s"*" SC state. The superconducting 
state is only weakly affected by doping. 

The goal of the present work is to understand how the 
system evolves from an SDW antiferromagnet to an s"*" 
superconductor. For this we derive and solve a set of 
coupled non-linear BCS-type equations for SC and SDW 
order parameters. We assume that the interactions in 
the two channels are comparable in strength and that 
Tc < r,, where Tg is the transition temperature at zero 
doping, Ts = Ts{x = 0). 

We report two results. First, when Tg/Tc is close to 
unity, the system displays second order SDW and SC 
transitions at Ts{x) and Tc, whichever is larger. The 
SDW state is commensurate, with momentum tt. At 
smaller T, the transition between SDW and SC upon 
changing x is first order, and there is no stable co- 
existence region (Fig. [T]). This is similar to the phase 
diagram reported for LaFeAsOi_a;Fa: in Ref[l|. Second, 
when Tg/Tc gets larger, SDW order becomes incom- 
mensurate with momentum Q = tt + q below some 
T* = 0.56r, > Tc (an SDW, phase-i^) We argue that 
in this situation SDWq and SC states co-exist. The co- 
existence region is initially confined to a small region 
below Tc, while at lower T the system still displays a 
first order transition between a commensurate SDW and 
SC states. As the ratio T^/Tc increases, the co-existence 
region extends down to lower T and eventually reaches 
T = (Figs. USD. 

An incommensurate SDWq state at finite dopings has 
been studied in connection with theoretical models for 
chromium and its alloys by Rica^ and othersji^ and in 
connection to pnictides by Cvetkovic and TcsanoviCfi^ 
Such state is a magnetic analog of Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) stated for which doping plays the 
same role as a magnetic field in a paramagnetically lim- 
ited superconductor.— We found that at T < Tc SDW^ 
phase exists only in combination with s+ superconduc- 
tivity. 
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Model and equations. We neglect double degener- 
acy of hole and electron states, which does not seem 
to be essential for the pnictides,— and consider a 
weak-coupling model with two families of fermions, near 
one hole and one electron FSs of small and near-equal 
sizes. The free electron part of the Hamiltonian is Ti.Q = 

('?c(k)4Q^Cka -hC/(k)/k„/kQ), where operators Cka 
and /ka describe fermions near (0,0) and (7r,7r), respec- 
tively (the momentum k in /ka is a deviation from tt). 



The dispersions ^■''''^(k) = ±^k+^, where ^k = vpik—kp), 
and S measures a deviation from a prefect nesting and 
may be tuned by doping(x) or pressure. The four-fermion 
part contains interactions in SDW and SC channels and 
in mean-field (BCS) approximation reduces to the ef- 
fective quadratic form H = 7 SkQ/3 ^kQ^^'k/3, with 

^'ka = (Cka, c_kQ, /k_|_qQ, /-k-qa), (^'- is a Conjugated 
column) : 
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The two diagonal blocks of the matrix H include the s+ 
SC order parameter A^ ~ — A-^ ~ A for two FS pockets 
and two off-diagonal blocks contain SDW parameter m^; 
^■^(k + q) = Ck + (5 + Vi?q for q <^ kp. The values of mq 
and A are determined by conventional self-consistency 
equations 



n^^V^'-'Y^al.ifl^^^c^p), (2a) 

k 

A = V'''Y.^-iny)^p{c.Mp) , (2b) 



where the sums are confined to only (0, 0) FS pocket, and 
ysdw ysc ^-^^ |j^g couplings in the particle-hole SDW 
channel and in the particle-particle SC s'^ channel.— 
Taken alone, V"^ leads to an s+ SC state with critical 
temperature Tc, while leads to an SDW state with 

transition temperature Ts at (5 = 0. The SDW order 
yields a real magnetization M{R) ~ cos QR at wave vec- 
tor Q = TT -I- q. The couplings and V^^ undergo 
logarithmic renormalizations from fermions with energies 
between ep and much larger bandwidth W and flow to 
the same value when W/ep 00 (Ref. ITsh . For any 
finite W/ep, V^"^"^ is the largest of the two. 

The correlators in Eqs. ^ are related to components of 
the Green's function G(k, r)a/3 = -(7V*(r)ka*(0)k/3), 



defined as the inverse of G 



lEn — where e„ 



7rT(2n -t- 1) are the Matsubara frequencies. The Green's 
functions in Eqs. ^ can be explicitly integrated over 
e- = [C^(k + q) - r(k)]/2 = ^k + ^v^q. Removing 
the coupling constants 2Nf\V'''\ and 2Nf\V'''^'"\ {Nf is 
the density of states at the Fermi surface per spin) and 
the upper cutoffs of the frequency sums in favor of the 



transition temperatures Tc^s, we obtain from Eqs. ([2]) 
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and 



In — = 2'kT ^ Re 




">0 

(3b) 

where angle brackets denote Fermi surface average. En ~ 
V4 + IAP, <5q = [^(k) + C/(k + q)]/2^ = 5 + 1 Vf^q, 
and Tg and T^, are solutions of the linearized equations, 
respectively for SDW (A = (5 = 0) and SC (m, = 0). 
This system of equations is solved to find all possible 
uniform SC and (generally) incommensurate SDW states. 

Note that for A = 0, Eq. (I3b| for SDW order coincides 
with the gap equation for a paramagnetically limited su- 
perconductor with niq instead of superconducting order 
parameter, 8 instead of magnetic field H , and incommen- 
surateness vector q instead of the total momentum of a 
Cooper pair ji^ii^ii^ 

We will also need the free energy F(A,mq) for these 
mean-field order parameters to pick out the state with 
the lowest F. We find F(A,TOq) in two complementary 
approaches: by differentiating with respect to interac- 
tion parameters^ and by using Luttinger-Ward func- 
tional and extending to a finite niq the derivation of the 
condensation energy for a BCS SCj ^°'^^ Both methods 
yield. 
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FIG. 1: (color online) (a) Phase diagram for Ts/Tc = 1.5 as 
a function of 5 controlled by doping. Thick solid and dashed 
lines are second-order SDW and SC transitions at Ta{5) and 
Tc, and dotted line - first order transition between commen- 
surate SDWo and s""" SC. Thin lines - physically unaccessible 
transitions. The pure magnetic Ts{5) line follows the curve 
of 'paramagnetically limited superconductivity'. The super- 
conducting Tc is independent of 5 in our model. Light lines 
denoted 5a and 5m are instability lines of SC and SDW states, 
(b) SDW and SC order parameters, and (c) free energy for 
SDW, SC and unstable mixed phases at T /Tc = 0.1. 



where AF(A,mq) = F{A,mq) - F(0,0). We solve self- 
consistency equations ([3]) for A and niq at finite S (i.e., 
doping) and arbitrary q and select the solution with q 
minimizing the free energy. 

Results. The results of our calculations are shown in 
Figs. [T]- [31 We find that the system behavior depends on 
the ratio Tg/Tc > 1. When this ratio is close to unity, the 
system only develops a commensurate SDW order with 
TOg=o = (see FigllJ. The SDW and SC transitions 



at Ts{S) and Tc (which is independent of 6, sec Eq. ([3a|) 
with mq = 0) are of second order. Below tricritical point 
at which Ts{S) = Tc, the transition between the states 
{m 7^ 0, A = 0) and (A 7^ 0, m = 0), is first order, and 
there is no region where m and A co-exist. 

The first-order transition at T = can be understood 
analytically. Setting q ~ and subtracting Eq. pap from 
(l3bD. we obtain for small 5: 



(5) 



Setting A^ = yields a linearized SC gap equations ([3]). 
We sec that in the presence of a nonzero mo = m(T ~ 



SI 



\HTslTc) 



Similarly, 

JC2 _ 



0), A first appears at 

for nonzero Aq, the SDW order grows from 
Ao ln(T's/rc). Their ratio, (Ja/^™ = too/Aq = Ts/Tc > 
1, implying that A nucleates in the SDW phase at a 
higher doping whereas m develops in SC state at lower 
6 (see FiglT})). This contradicts the very fact that SDW 
state is stable at smaller dopings than the SC state. The 
solution with A, m 7^ then grows in "wrong" direction 



FIG. 2; (color online) Same as in Fig. [T] but for Ts/Tc = 5. 
An incommensurate SDW order appears below Ts{S) once 
it becomes smaller than T* — 0.56Ts > Tc- Below Tc, a 
new mixed phase appears in which incommensurate SDW, 
order co-exists with SC. At small T, there is no SDW, state 
without superconductivity. The transitions into the mixed 
state are second order from a SC state and first order from a 
commensurate SDW state. The free energy now shows that 
near 5/2ttTs = 0.2 a mixed state has lower energy than the 
two pure states. 



of S, and we explicitly verified that is has a higher energy 
than pure states and therefore is unstable, see Fig. [TJ:. 
As both q = SDW and SC gaps cover the entire FS, 
the absence of the state where the two co-exist implies 
that fully gapped SC and SDW orders cannot co-exist, 
and only one of these two states is present at a given S. 
First order transition between the SDW and SC states 
occurs at 6 ^ Scr, when their free energies coincide. This 



happens when —711^/2 + 5^ 



-A^/2, hence 



A2 = 

2 



1- ;^ 



(6) 



which is in between the two second order instability 
points 5m and 5/^. 

Situation changes when Tg/Tc gets larger, and there 
appears a wider range of dopings where Ts{5) > Tc- If 
only commensurate magnetic order SDWq was possible, 
magnetic transition would become first-order below a cer- 
tain T* = O.SGTs, which at large enough Ts/Tc becomes 
greater than Tc. In reality, the system avoids a first-order 
transition and extends the region of magnetic order by 
forming an incommensurate SDWg state below T* (see 
Figs. [21 [3]) fi^ii^ii^ The transition from the normal state 
to SDWg state is second order, the subsequent transition 
to the commensurate SDWq state is first order. Once nig 
is developed, the actual solution is more complex and in- 
cludes harmonics with multiple q (Ref. [l^ . but we ignore 
this for now. 

Our main result is the discovery of a new phase below 
Tc, in which incommensurate SDWg order co-exists with 
SC order. Physically, the key reason for appearance of 
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FIG. 3: (color online) Same as in previous two figures, but for 
intermediate Ts/Tc = 3. Mixed phase appears only in a tiny 
region near the point where Tc and Ts {5, q) cross. At smaller 
T, the system still displays a first order transition (dotted 
line) between a commensurate SDWo state and a SC state. 



such phase is that incommensurate SDWg order does not 
gap the excitations on entire FS ~ the system remains a 
metal albeit with a modified, smaller FS.— i^i^ Once the 
Fermi surface survives, an attractive pairing interaction 
gives rise to SC below Tc. Alternatively speaking, for 
SDWg order, some parts of the FS become unaccessible 
to magnetic 'pairing', and SC order takes advantage of 
this, c.f. Refs. mil 

At large enough Tg/Tc, the co-existence phase extends 
to T = (see Fig. [2] for T^/Tc = 5). In Fig. Ha we show 
order parameters at T = Tc/2. The transition from a 
commensurate SDW state into a mixed state is first order 
with both A and the amount of incommensurateness q 
jumping to finite values. The transition from a SDWq 
state into a mixed state, as well as the transition from a 
mixed state into a pure SC state are of second order with 
ruq gradually vanishing. Fig.[2t shows the corresponding 
free energies of all states. Comparing it with Fig. [T]; we 
see that now the SC state becomes unstable at a higher 
5 and the mixed state now has lower energy than pure 
SC or SDW states. 

The behavior at somewhat smaller Ts/Tc is intermedi- 
ate between those in Figs. [l]and [H In figure [3] we show 
the phase diagram for intermediate Tg/Tc = 3. We still 
have T* > Tc, and the mixed phase still exists, but it 



now appears only as a small pocket near the point where 
Ts{S) = Tc- At smaller T, the system shows the same 
behavior as at q = 0, i.e., a first-order transition between 
commensurate SDW and SC states. 

To conclude, in this paper we considered SDW and s+ 
SC instabilities in a mean-field approximation for a two- 
band model for the pnictides. We assumed that there 
are "attractive" interactions in both SDW and extended 
s-wave SC channels, leading to the nesting-driven SDW 
magnetism at Ts{S), decreasing with S controlled by dop- 
ing or pressure, and to s+ SC at T^. At zero doping, 
Ts(0) ~ Ts > Tc, but Tc becomes the first instability at 
large dopings. The issue we consider is how SDW or- 
dered state transforms into a superconducting s+ state. 
We found that the transition is first order, with no in- 
termediate mixed phase, if Tg and Tc are close enough 
so that SDW order above Tc is commensurate. At larger 
Ts/Tc, the system develops an incommensurate SDWg or- 
der above Tc- This incommensurate SDW order does not 
gap the whole FS and allows for a co-cxistcncc of mag- 
netism and superconductivity. We found that the mixed 
phase first appears in a small pocket near Ts{6) = Tc, but 
extends as Ts/Tc increases eventually reaching down to 
T = 0. At even larger Tg/Tc, zero temperature phases in- 
clude a commensurate SDW phase, an s'^ SC phase, and 
an intermediate mixed phase where SC co-exists with 
an incommensurate SDW order, but there is no purely 
incommensurate SDW at T = 0. The transition from 
the mixed phase to a commensurate SDW phase is first 
order, and to the SC phase is second order. 

Finally, in this paper, we only considered a mixed state 
with a uniform SC order (zero total momentum of a pair). 
In principle, an inhomogencous SC state (a true FFLO 
state) is possible because incommensurate SDW^ order 
breaks the symmetry between FS points with k and — k. 
This, however, should not change the phase diagram as 
a non-uniform SC state may only exist at large enough 
•niq, i.e., near a first order transition into a commensurate 
SDW phase, possibly extending a mixed state into SDWq 
region. For large rUq, an approximation by a single q 
is not sufficient. A more sophisticated analysis in this 
region is called for. 
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